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The proton-proton momentum correlation function (Cpp) from different rapidity regions are sys- 
tematically investigated for the Au + Au collisions at different impact parameters and different 
energies from 400A MeV to 1500A MeV in the framework of the isospin-dependent quantum molec- 
ular dynamics model (IQMD) complemented by the Lednicky and Lyuboshitz analytical method. 
In particular, in-medium nucleon-nucleon cross section (NNCS) dependence of the correlation func- 
tion is brought into focus, while the impact parameter and energy dependence of the momentum 
correlation function are also explored. The sizes of the emission source are extracted by fitting the 
momentum correlation functions using the Gaussian source method. We find that the in-medium 
nucleon-nucleon cross section obviously influence the proton-proton momentum correlation function 
which is from the whole rapidity or projectile/target rapidity region at smaller impact parameters, 
but there is no effect on the mid-rapidity proton-proton momentum correlation function, which 
indicates that the emission mechanism differs between projectile/target rapidity and mid-rapidity 
protons. 


PACS numbers: 25.70.Mn, 24.10.-i, 25.70.Pq, 27.80.--w 


I INTRODUCTION sides the applications of the HBT method to investigate 


The Hanbury Brown and Twiss (HBT) effect was first- 
ly discussed in radio and stellar astronomy. The method 
was applied to measure the angular diameter and the size 
of stars by Hanbury Brown and Twiss [1]. Later on, the 
technique was introduced to research particle physics in 
1960s by Goldhaber et al.. They studied the angular 
distribution of identical pion pairs in proton-antiproton 
annihilations and observed an enhancement of pairs at 
small relative momenta [2]. In the last decade, there has 
brought great strides in experiments and a larger number 
of theoretical researches ranging from low energy to high 
energy HICs [3, 4]. It is well known that the two-particle 
correlation is sensitive to characteristics of the particle 
emission source. Recently, the two particle correlation 
in subatomic physics has been taken as a probe for the 
space-time geometry of the particle emission source. Cor- 
relation between two protons was measured by several 
experiments and explored by different models. Not only 
protons, but also the composite light fragments/particles, 
which will not be discussed in this paper, are also used to 
carry information on the emission source [5, 6]. Due to 
the rapid development of radioactive nuclear beams, the 
HBT method is also used to study the exotic structure 
of nuclei. For instance, there have been several measure- 
ments for revealing exotic structure of the neutron-rich 
nuclei such as 9He, !! Li, !4Be [7-9] and of proton-rich nu- 
clei such as ?? A] [10] as well as ??Mg [11, 12]. In addition, 
the dependence of the proton-neutron correlation on the 
binding energy was also theoretically explored [13]. Be- 


* Corresponding author: ygmaQsinap.ac.cn 


the exotic structure, it has become an important tool in 
heavy-ion collisions at wide energy range [3, 4, 14-16]. 
For example, in relativistic energy region, the Collabora- 
tions at the RHIC and the LHC have carried out a lot 
of experimental measurements of the correlation func- 
tion of two-pion as a function of energy and system size 
[17, 18]. What's more, the same method has been ap- 
plied to make the first measurement of two-antiproton 
interaction by analyzing momentum correlation function 
between antiprotons, namely the quantitative extraction 
of the scattering length and the effective range which are 
two key parameters to characterize the strong interaction 
for the antiproton interactions by the STAR Collabora- 
tion [19-22]. Theoretically, the correlation functions be- 
tween two identical pions or kaons were also investigated 
in some simulation work such as the hydrodynamic mod- 
el and the AMPT model etc [4, 23-27]. In intermediate 
energy region, the two-proton correlation functions have 
been mostly applied to extract the space-time properties 
such as the source size and emission time in the nuclear 
reaction [28]. In addition, there are many investigation- 
s on the dependences of the correlation functions in the 
experiments and theories, such as on impact parameter 
[29, 30], total momentum of nucleon pairs [31], isospin of 
the emission source [32], nuclear symmetry energy [33], 
nuclear equation of state (EOS) [30], the density distri- 
bution of valence neutrons in neutron-rich nuclei [6] and 
so on. 


The investigation of in-medium nucleon-nucleon scat- 
tering is of interest in intermediate energy heavy ion re- 
actions. In this energy domain, nucleus-nucleus collisions 
provide a unique opportunity to form the compressed nu- 
clear matter with a density up to 2-3 times normal nucle- 
ar matter density (9). The in-medium nucleon-nucleon 


cross section (NNCS) has a close relation with the nuclear 
matter density. Therefore, it is an important component 
in our model simulations. Recently, the medium effects 
on nucleon-nucleon cross section have been widely inves- 
tigated by replacing the NNCS in vacuum with an in- 
medium one and the various effects have been discussed 
[35-38]. 

'The dependence of the two-proton correlation function 
on the in-medium NNCS has been briefly studied through 
CRAB code in a framework of IQMD [30]. Since two- 
particle correlation function, through final-state interac- 
tions and quantum statistical effects, have been shown 
to be a sensitive probe to the space-time distributions of 
emitted particle in heavy-ion collisions [16], it is of great 
interest to investigate the in-medium NN cross section 
effects on the source evolution. In the present paper, we 
use another theoretical approach which was proposed by 
Lednicky and Lyuboshitz [39] to explore the relationship 
between the above factors and proton-proton correlation 
function more details. The two-particle correlation at s- 
mall relative velocities is sensitive to the space-time char- 
acteristics of the production process owning to the effects 
of quantum statistics and final-state interaction [40, 41]. 
In most proton-proton correlation functions, the HBT 
strength at 20 MeV/c of the p-p relative momentum is 
taken as a unique quantity to determine the source size 
or emission time of two-proton emission [30]. The proton 
phase spaces of Au+Au collisions at the freeze-out time 
generated by the IQMD model are used as the input for 
the Lednicky and Lyuboshitz code and then the effective 
source size of the source will be extracted. 

The rest of paper is organized as follows. In Sec. II 
we briefly describe the models and formalism used in the 
present study, i.e., Lednický and Lyuboshitz analytical 
formalism and an isospin-dependent quantum molecular 
dynamics model (IQMD) model. The detailed analysis 
and discussion of systematic proton-proton momentum 
correlation function and extracted source size results for 
different rapidity regions are given in different in-medium 
nucleon-nucleon cross sections, different impact parame- 
ters, beam energies for Au + Au collisions at in Sec. III. 
In addition, we fit the proton pr spectra with the distri- 
bution function from the Blast-Wave model. Finally, in 
Section IV we summarize the results. 


IIl. FORMALISM AND MODELS 


A. LEDNICKY and LYUBOSHITZ 
ANALYTICAL FORMALISM 


Firstly, we would like to present a brief review of the 
theoretical approach which was proposed by Lednicky 
and Lyuboshitz [39] for the HBT analysis. The method 
is based on the principle that correlation function of i- 
dentical particles when they are emitted at small rela- 
tive momenta are determined by the effects of quantum- 
statistical symmetry (QS) of particles and the final- 


state interaction (FSI) [42]. In this technique, we as- 
sume the particles emitted by independent one-particle 
point source and the spin independent in the production 
progress as well as the two-particle interaction. Then we 
can investigate particle pairs (1,2) emitted at small rela- 
tive momenta. Due to the conditions following references 
[43] , we neglect the effect of FSI in all pairs (1,2) and 
(2, 1) except (1,2). We can see the progress in fig 1. 
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FIG. 1: The diagram describes the production of interacting 


particle 1 and 2 [43]. 


So the correlation function for identical particles takes 
the expression: 


B (p,q) = Bo (p,q) + Bi (p,q), (1) 
where 
1 
P=Pi+P2,4= 5|P2— Pil (2) 


are the total momentum and relative momentum of the 
particle pair, respectively. In Eq. (1), Bo (p,q) is the 
contribution of quantum statistics effect, described by 
the formula: 


Bo (p. 4) = go cos (qx) . (3) 


Here go is the spin factor. 

Then, the function B, (p,q) can be expressed through 
the symmetrizied Bethe-Salpeter amplitude 4% (S) which 
can be approximated by the outer region solution of the 
scattering problem [19]: 


Bı (p,q) = S psp (21,22, 23, 25) x VS a (21,22) 
S 


x VS a (21,25) d'zid'zad'zd!z), (4) 


with psp is the two-particle density matrix. 

Next, we can introduce the detailed analytical calcula- 
tion of the proton-proton correlation function [19]. The 
proton-proton correlation function, Cpp (k*, ro) can be 
described by the Lednicky and Lyuboshitz analytical 
method [19, 39]. In this model, the space distribution 
of the Gaussian source was simulated according to the 
following function: 


S (r*) & exp (-r""/ (4r3) ) i (5) 


Here, ro is the source size parameter. Therefore, we can 
obtain the correlation function through assuming 1/4 of 
the singlet and 3/4 of triplet states. The theoretical cor- 
relation function at a given k* can be calculated as the 
average FSI weight (w(k*,r*)) obtained from the sepa- 
ration r*, simulated according to the Gaussian law, and 
the angle between the vectors k* and r*, simulated ac- 
cording to a uniform cosine distribution. The average 
FSI weight can be described by the formula: 

SCH) c." 

e (t")) /2 (6) 


w (k*,r*) = 9? C? (r* 


)+(-1)° 
where S is the total pair spin, r* is the relative distance, 
PM (r*) is the equal-time (t* — 0) reduced Bethe- 


Salpeter amplitude which can be approximated by the 
outer solution of the scattering problem [39]. This is 


Rr (r*) = eJ A. (A)X 


peus (—iA, 1, i£) + fe (k*) CoN) x. A 


where ôe —argI (1 + i/k*a;) is the Coulomb phase cor- 
responding to zero orbital angular momentum, A, (A) = 
2A [exp (21A) — 1]! determines the contribution of the 
Coulomb interaction, i.e., the positive value correspond- 
ing to the repulsion, \ = (k*az) t, ac = 57.5 fm is the 
Bohr radius for two protons, p = k*r*,£ = k*r* + p, 
F is the confluent hypergeometric function, G (p, A) = 
VA (X) [Go (p, A) + iFo (p, À)| is a combination of the 


Ed i and singular (Go) s-wave Coulomb function, 


—1 


ik'A,Q)! (8) 


fane jd" - —h(d) 
is the s-wave scattering amplitude renormalizied by the 
Coulomb interaction, do is the effective range of the inter- 
action, and A (A) = A? OP, [n (n? + PS —C —ln[A 
(here C = 0.5772 is the Euler constant). The depen- 
dence of the scattering parameters on the total pair spin 
S is omitted since only the singlet (S — 0) s-wave FSI 
contributes in the case of identical nucleons. 


B. THE IQMD MODEL 


To apply the above theoretical simulation, the single- 
particle phase-space distribution at the freeze-out is re- 
quired. In this work, the correlation function can be 
established from the emission phase space given by the 
IQMD transport model [44] . 

The quantum molecular dynamic model (QMD) is a 
many-body transport theory, it has been extensively ap- 
plied to describe heavy-ion reactions from intermediate 
energy to 24 GeV [45]. By the QMD studies, various 
valuable information about both the collision dynamics 
and the fragmentation process has been learned [46-55]. 


Excellent extensibility can be also expected due to its mi- 
croscopic treatment on the collision process. The model 
mainly consists of several parts: initialization of the pro- 
jectile and the target nucleons, nucleon transport under 
the effective potentials, nucleon-nucleon (NN) binary col- 
lisions in a nuclear medium, the Pauli blocking, and the 
numerical test. The isospin-dependent quantum dynamic 
model (IQMD) is based on the QMD model and consid- 
ers the isospin factors [56] in mean-field, two-body NN 
collisions, and the Pauli blocking. In the IQMD mod- 
el, the wave function of each nucleon is represented by 
the form of Gaussian wave packet, with the parameter L 
which relates to the size of the reaction system. For Au 4- 
Au system, the width L is fixed to 2.16 fm? . Gaussian 
wave packet is written as 


ói(r) = 
1 -(r-ri() ir - p, (t) 
a2nL5/4 ^P | AL nd ( h ) e 


Here, r; (t) and p; (t) are the time dependent variables 
which describe the center of the packet in coordinate and 
momentum space, respectively. Then, all nucleons inter- 
act via the effective mean-field and two body NN colli- 
sions. 

The nuclear mean field can be expressed as: 


U = Usky +UcCoui t+Uvur+Usymt+Umpi+UPaui, (10) 


where Usky, UCoul, UyYuk, Usym, Umpr, and Upaui are 
the density-dependent Skyrme potential, the Coulom- 
b potential, the surface Yukawa potential, the isospin 
asymmetry potential, the momentum-dependent interac- 
tion and the Pauli potential, respectively. A general re- 
view of the above potentials can be found in Ref. [45]. 
In the present work, the in-medium NN cross section is 
represented by the formula: 


oped = (1 : ae) ofr, (11) 


where po is the normal nuclear matter density, p is the 
local density, r is the in-medium factor and of"<* is the 
available experimental NN cross section [57]. In this 
above expression, increasing values of parameter 7 cor- 
respond to decreasing values of the in-medium nucleon- 
nucleon cross-section. 

In this model, the fragments are identified using a mod- 
ified minimum spanning tree description. In the Mini- 
mum Spanning Tree approach, two nucleons are assumed 
to share the same cluster if their centers are closer than 
a distance of 3.5 fm and their relative momentum small- 
er than 0.3 GeV/c. If nucleon is not bounded by any 
clusters, it is treated by an emitted (free) nucleon. 


C. THE BLAST-WAVE FIT 


In heavy-ion collisions, particles collide with each other 
randomly, which can be described in term of thermal mo- 


tion [58]. We adopt the blast-wave model which has been 
put forward by Siemens and Rasmussen [59] to describe 
the mid-rapidity pr spectra with two free parameters: 
collective transverse flow velocity 8 and kinetic freeze-out 
temperature Tf. The collective transverse flow velocity 
D is parametrized by the surface velocity 6, in the region 
of 0 € R< Ras: 


8i) - (qr) . (12) 


where Rmax is the maximum radius of the expanding 
source at thermal freeze-out time, and the 6, is the 
particle radial velocity at the maximum surface, e.g. 
r = Rmazx, and the exponent o describes the evolu- 
tion of the flow velocity with the radius. The pr spec- 
tra are a superposition of individual thermal sources 
with different r, which is boosted with the boost angle 
p = tanh”! B, (r) [60? ]. 
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where, Kı, Io are the modified Bessel functions. The 
shapes of the spectra are essentially determined by TY, 
Bs, a, and the mass of the particle mo. The average 
flow velocity is estimated by taking an average over the 
transverse geometry. 


Ill. ANALYSIS AND DISCUSSION 


In this work, we use the soft EOS with momentum 
dependent interaction for all Au + Au collisions at beam 
energies from 0.4 to 1.5A GeV. The correlation functions 
are calculated by using the phase-space information from 
the freeze-out stage. 

Firstly, we investigate the influence of the in-medium 
NN cross section on the momentum correlation function 
for Au + Au collisions at 1.04 GeV. Fig. 2 shows the 
proton-proton momentum correlation function for Au + 
Au collisions at 14 GeV with different in-medium re- 
duction factors 7, impact parameters, and proton rapid- 
ity region. In each panel, in-medium reduction factor 
of 0.0, 0.2, 0.5, 0.7 and 0.9 are compared. From top 
panel to bottom panel, each panel corresponds to dif- 
ferent impact parameters, from b = 3, 6, 9 up to 12 
fm, respectively. From left to right columns, it repre- 
sents the correlation function for the proton pairs with- 
in whole rapidity, mid-rapidity and projectile/target ra- 
pidities, respectively. Here the mid-rapidity cut means 
that both protons are emitted in the rapidity window of 
-0.5< y/yproj < 0.5, and projectile or target rapidity 
region means that both protons either come from the ra- 
pidity region of y/Yproj. > 0.5 or y/yproj < —0.5, where y 


represents the proton rapidity and Yproj. means the ini- 
tial projectile rapidity. Overall, the proton-proton mo- 
mentum correlation function exhibits a peak at relative 
momentum q = 20 MeV/c, which is due to the strong 
final-state s-wave attraction together with the suppres- 
sion at lower relative momentum as a result of Coulomb 
repulsion and antisymmetrization wave-function between 
two protons. 


For protons which are emitted in whole rapidity or pro- 
jectile/target rapidity, the general trend is very similar. 
With the increasing of the in-medium NN cross section 
(i.e. the less in-medium reduction factor 7), the collision 
rate between nucleons increases. Therefore more nucle- 
ons are emitted early, which make the strength of the 
momentum correlation function larger. The difference 
can be further revealed in central and semi-peripheral 
collisions, however, the difference in momentum correla- 
tion function among different r-factors almost disappears 
in peripheral collisions. This indicates that the NN cross 
section in peripheral collisions has no significant change 
even though the 7 value changes much. Overall speaking, 
with the less nucleon-nucleon cross section (i.e. larger 7 
factor), the correlation peak decreases, indicating that 
proton-proton correlation has a positive correlation with 
nucleon-nucleon cross section. The sensitivity of corre- 
lation strength to the r values becomes less important 
when the reaction goes to peripheral collisions. For mid- 
rapidity protons, correlation functions are much stronger 
than the cases of whole-rapidity or project/target ra- 
pidity, and shows almost no dependence on in-medium 
nucleon-nucleon cross section, which indicates very dif- 
ferent spacial-time structure of mid-rapidity protons. Es- 
sentially, the mid-rapidity protons are emitted very ear- 
ly and non-equilibrium, therefore insensitive to the in- 
medium NN cross section, additionally they have radial 
flow which will be mentioned later. In addition, another 
difference is the strength of correlation peak as a function 
of impact parameter is reverse to the behavior of correla- 
tion peak for whole rapidity or projectile/target rapidity 
protons, which displays a stronger correlation peak in 
peripheral collision than central collisions, indicating s- 
lighter more compact mid-rapidity source in peripheral 
collisions. 

How the three variables including the incident energy, 
impact parameter and in-medium nucleon-nucleon cross 
section affect the strength of Cpp is presented in Fig. 3 
for proton-proton correlation function where the whole 
rapidity window (blue squares) or mid-rapidity (solid cir- 
cles) is applied for emitted protons in Au + Au colli- 
sions. For the incident energy dependence, the increase 
of momentum correlation peak from low energy to high 
energy is cleanly seen in panel (a) of Fig. 3 , and the 
mid-rapidity cut displays a stronger peak. It can be gen- 
erally understood that a more rapid collision process, and 
a smaller emission source space and short time interval 
among emitting nucleons at higher energy, so does the 
mid-rapidity source [30]. In the panel (b) of Fig. 3 , we 
see that the correlation strength increases with the im- 


pact parameter for the mid-rapidity cut but opposite for 
the whole rapidity cut, which shows that geometrical cut 
for two rapidities are complemented, e.g. stronger cor- 
relation for mid-rapidity protons at peripheral collision 
indicates smaller source size, but weaker correlation for 
whole-rapidity protons indicates larger source size. In the 
panel (c) of Fig. 3 , with the increase of the in-medium 
cross section modification factor, i.e. decreasing the in- 
medium NN cross section, the peak strength becomes 
smaller. In the other words, proton-proton correlation 
positively depends on the nucleon-nucleon collisions. On 
the other hand, in contrast with the whole-rapidity pro- 
tons, mid-rapidity protons show weaker sensitivity to the 
NNCS. As we have already known that the strength of 
the correlation function depends mainly on the source 
size, the above behavior of the HBT strength essential- 
ly reflect the changing size of the emission source versus 
the beam energy, impact parameter and nucleon-nucleon 
Cross section. 


Fig. 4 presents the radius of the Gaussian source for 
emitted protons as a function of beam energy (left colum- 
n), impact parameter (middle column) and in-medium 
cross section reduction factor (right column) for differ- 
ent rapidity windows, namely the whole rapidity (upper 
row), mid-rapidity (middle row) and projectile/target ra- 
pidity (bottom row). Overall, for mid-rapidity proton- 
proton correlations, the sensitivity to the in-medium 
cross section reduction factor is almost not visible. The 
source radius shows a slight drop with the the increasing 
of beam energy or impact parameter. For the whole ra- 
pidity or projectile/target rapidity windows, their beam 
energy dependences are very similar if the same 77 is ap- 
plied, and source radius increases with the increasing of 
1, ie. the decreasing of nucleon-nucleon cross section. 
When the nucleon-nucleon cross section is larger (e.g. 
1] € 0.5), source radius drops with the incident energy, 
indicating fast emission and/or compact proton emission 
size in higher energies. However, the situation is differ- 
ent when the nucleon-nucleon cross section is small (e.g. 
n > 0.7), where radius shows a slight increasing or a 
plateau behaviour. But overall, even though slight beam 
energy dependence of radius is seen, the trend is rather 
weak as seen in ultra-relativistic energy heavy ion colli- 
sions [63]. 

Middle column and right column demonstrate the ra- 
dius of the Gaussian source as a function of impact pa- 
rameter at different fixed 7 values or 7 at different im- 
pact parameters, respectively. As expected, for corre- 
lation between protons from the whole rapidity region 
or projectile/target region, the source size increases with 
the increasing of impact parameter, and the larger the 
nucleon-nucleon cross section, the stronger the depen- 
dence of the source size on impact parameter. With 
increasing impact parameters, the effective source size 
gets bigger because apparently the (target or projectile) 
spectator region gets bigger and bigger with increasing 
impact parameter, and this effect is more pronounced 
if the in-medium decrease of the cross-section is small- 


er. This indicates some geometrical effect and protons 
coming from spectator fragmentation mechanism. At the 
same time, the source size increases with the decreasing 
of the in-medium NN cross section (i.e. larger 7 factor). 
'The smaller the impact parameter, the stronger the de- 
pendence of the source size on in-medium NN cross sec- 
tion, which can be understood that more frequent colli- 
sion effect is expected to affect the dynamical evolution. 
However, for correlation between protons from the mid- 
rapidity region, the tendency of the source size becomes 
decreasing with the increasing of impact parameter, re- 
gardless of the in-medium nucleon-nucleon cross section. 
'This reverse dependence of source radius between mid- 
rapidity protons and projectile/target rapidity protons as 
a function of impact parameter indicates some geometri- 
cal evolution of participant and spectator region, which 
can be also well seen from Fig. 3(b). 


We should mention that to extract the above source 
size, theoretical calculations for Cpp was performed by 
using the Lednicky and Lyuboshitz analytical method. 
The best fitting source size is judged by finding the min- 
imum of the reduced chi-square. Fig. 5 present examples 
of the y?-variance between the IQMD calculations with 
the Lednicky-Lyuboshitz analytical formalism and the 
Gaussian source correlation as a function of the radius of 
the Gaussian source in different impact parameters with 
the changing 7 factors (left panels) or in different 7 fac- 
tors with the changing impact parameters (right panels). 
Generally, the minimum can be well defined, but the er- 
rors on its location are apparently rather asymmetric. 


Middle rapidity protons are possibly experienced by 
the collective radial flow expansion in comparison with 
the projectile/target rapidity protons. To demonstrate 
how large the radial flow is for the mid-rapidity protons, 
we use the Blast-wave (BW) fits to pr spectra of protons 
in the mid-rapidity region. Here, we get the a value on 
the pr spectra of central collisions and fix this value when 
do the fitting on the other impact parameters [34]. Left 
column of Fig 6 shows the pr spectra of mid-rapidity pro- 
tons in different incident energy for b = 3 fm and ņ = 0.2 
(top row), in different impact parameter for E = 1A GeV 
and 7 = 0.2 (middle row) and with different in-medium 
cross section for b = 3 fm and 1A GeV (bottom row), 
where the solid symbols represent the calculated results 
from the IQMD and the solid lines are the BW fits. Over- 
all all lines can well reproduce the spectra, from which 
the radial flow parameters can be systematically extract- 
ed. Right column of Fig. 6 displays the extracted radial 
flow parameters (8) as a function of incident energy at b 
= 3 fm and 7 = 0.2 (top panel), as a function of impact 
parameter at E = 1A GeV and 7 = 0.2 (middle panel), 
and as a function of in-medium cross section factor at E 
— 1A GeV and b — 3 fm (bottom panel). Obviously, the 
radial flow becomes stronger in higher incident energy 
as well as in more central collisions. Meanwhile, larg- 
er in-medium nucleon-nucleon cross-section (i.e. smaller 
7 values) leads to larger radial flow due to the frequen- 
t nucleon-nucleon collision in overlap zone. It is noticed 


that radial flow has been already extensively discussed in 
mediate and high energy HIC, eg. discussed in Ref. [62] 
by J.Helgesson et al.. 

For mid-rapidity p-p correlation, let us have a close 
look for correlation strength versus the radial flow. Since 
we have relationship between the correlation strength at 
20 MeV/c (Cpp(q=20MeV/c) ) versus beam energy as 
well as the radial flow (8) versus beam energy, we can 
obtain the relationship between C, (q—20MeV /c) and 8 
which is displayed in Fig. 7(a) in the condition of b — 
3 fm and 7 = 0.2. It tells us that the larger the radial 
flow velocity, the stronger the proton-proton correlation. 
In the same way, we got Cpp(q=20MeV/c) versus B in 
Fig. 7(b) in the condition of E = 1A GeV and 7 = 0.2 
where the impact parameter is a variable. In this case, 
anti-correlation is observed. Simiarily, C5, (q220MeV /c) 
versus f is shown in Fig. 7(c) in the condition of E=1A 
GeV and b = 3 fm where the in-medium nucleon-nucleon 
cross section reduction factor is a variable. Here, a s- 
light increasing behavior is demonstrated. Seen from the 
above three variables, we found no unique dependence 
of p-p correlation function as a function of radial flow 
parameter. 


IV. SUMMARY 


In the present work, we use the IQMD transport ap- 
proach to calculate the phase-space points at the freeze- 
out stage for Au + Au collisions from 0.4 to 1.5A GeV. 
Afterwards the phase space points were processed within 
the Lednicky and Lyuboshitz analytical formalism to re- 
construct the proton-proton correlation function. In this 
way, we systematically study how the in-medium NN 
cross section affects the strength of the momentum cor- 
relation function of proton-proton pairs for Au + Au from 
0.4 to 1.5A GeV in different rapidity windows and dif- 
ferent impact parameters. Results show that the larger 
in-medium NWN cross section results in a stronger mo- 
mentum correlation function than a smaller in-medium 
NN cross section, especially at small impact parame- 
ters, for the whole rapidity or projectile/target rapidity 
proton pairs. This behavior is interpreted as stronger 
correlation for equilibrium-like protons induced by the 
higher nucleon-nucleon collision rate. However, for the 
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mid-rapidity proton emission, the in-medium nucleon- 
nucleon cross section has less effect on the momentum 
correlation function due to very different emission mech- 
anism. In addition, the impact parameter effect on the 
HBT strength has been also addressed in the work. We 
have shown that the HBT strength has very different de- 
pendence between the whole rapidity or projectile/target 
rapidity proton pairs and the mid-rapidity proton pairs. 
For proton pairs from projectile/target rapidity, it de- 
creases with the increasing of impact parameter; but for 
proton pairs from mid-rapidity, it increases with the in- 
creasing of impact parameter. By fitting the momentum 
correlation function with the Gaussian source, the effec- 
tive proton emission source sizes are extracted. Results 
show that the source radius general increases with the 
increasing of the impact parameter for emitted protons 
within the whole rapidity or projectile (target) rapidi- 
ty, however, it decreases versus the impact parameter 
for the mid-rapidity protons and it shows insensitivity to 
the in-medium nucleon-nucleon cross section. The above 
phenomenon reflects the evolution of source size with the 
collision geometry, i.e. the mid-rapidity source becomes 
smaller but the projectile/target source becomes larger 
with the increasing of impact parameter. Moreover, the 
beam energy dependence of HBT strength is also present- 
ed. Generally, the HBT strength shows a slight change 
with beam energy especially for the whole rapidity or 
projectile (target) rapidity proton pairs. By using the 
Blast-wave fits to the transverse momentum spectra of 
mid-rapidity protons, the radial flow parameters are sys- 
temically extracted as a function of beam energy, impact 
parameter and in-medium nucleon-nucleon reduction fac- 
tor, and therefore relationships of HBT strength versus 
radial flow parameter are constructed in different condi- 
tions. However, no unique dependence is found, which 
indicates the radial flow is not a decisive variable for the 
p-p correlation. 


This work was supported partly by the National 
Natural Science Foundation of China under Contract 
Nos. 11421505, 11220101005, the Major State Basic Re- 
search Development Program in China under Contract 
2014CB845401, and the Key Research Program of Fron- 
tier Sciences of the CAS under Grant No. QYZDJ-SSW- 
SLH002. 


6| X. G. Cao, Y. G. Ma, D. Q. Fang et al., Phys. Rev. C 
86, 044620 (2012). 

7] K. Ieki et al., Phys. Rev. Lett.70, 730 (1993); F. M. 
Marques et al., Phys. Lett. B 476, 219 (2000). 

8| M. T. Yamashita, T. Frederico, and L. Tomio, Phys. Rev. 
C 72, 011601(R) (2005); M. Petrascu et al., Phys. Rev. 
C 69, 011602(R) (2004). 

9] Z. Kohley et al., Phys. Rev. Lett. 110, 152501 (2013); E. 
Lunderberg et al., Phys. Rev. Lett .108, 142503 (2012). 
[10] P. Zhou et al., Int. J. Mod. Phys. E 19, 957 (2010); X. 


11 


12 


13 


14 


15 


16 


17 


18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 


34 


Y. Sun et al., Int. J. Mod. Phys. E 19, 1823 (2010). 

Y. G. Ma, D. Q. Fang, X. Y. Sun et al., Phys. Lett. B 
743, 306 (2015). 

D. Q. Fang, Y. G. Ma, X. Y. Sun et al., Phys. Rev. C 
94, 044621 (2016). 

Y. B. Wei et al., Phys. Lett. B 586, 225 (2004); Y. B. 
Wei et al., J. Phys. G 30, 2019 (2004). 

S. Pratt and M. B. Tsang, Phys. Rev. C. 36, 2390 (1987); 
S. E. Koonin, Phys. Lett. B 70, 43 (1997); S. Pratt, Phys. 
Rev. Lett. 53, 1219(1984). 

J. P. Sullivan et al., Phys. Rev. Lett. 70, 3000 (1993). 
W. Bauer, C.K. Gelbke, and S. Pratt, Annu. Rev. Nucl. 
Part. Sci. 42, 77 (1992). 

C. Adler et al. (STAR collaboration), Phys. Rev. Lett. 
87, 082301(2001); J. Adams et al. (STAR collaboration), 
Phys. Rev. Lett. 93, 012301 (2004); J. Adams et al. (S- 
TAR collaboration), Phys. Rev. Lett. 71, 044906(2005); 
S. S. Adler et al. (PHENIX collaboration), Phys. Rev. 
Lett. 93, 152302 (2004); B. Abelev et al. (STAR collabo- 
ration), Phys. Rev. C 80, 024905 (2009); B. Abelev et al. 
(STAR collaboration), Phys. Rev. C 81, 024911 (2010). 
K. Aamodt et al. (ALICE collaboration), Phys. Lett. B 
696, 328 (2011). 

L. Adamczyk et al. (STAR Collaboration), Nature 527, 
345 (2015). 

Z. Q. Zhang, Ph.D. Dissertation, Shanghai Institute of 
Applied Physics, Chinese Academy of Sciences, China, 
May 2017. 

H. Q. Zhang, Nat. Sci. Rev. 3, 154 (2016). 

Z. Q. Zhang, Y. G. Ma, Nucl. Sci. Tech. 27, 152 (2016). 
Zheng-Qiao Zhang, Song Zhang, Yu-Gang Ma, Chin. 
Phys. C 38, 014102 (2014). 

Z. W. Lin et al., Phys. Rev. Lett., 89, 152301 (2002); J. 
Phys. G: Nucl. Part. Phys. 30, $263 (2004). 

F. Retiére and M. Lisa, Phys. Rev. C 70, 044907 (2004). 
M. A. Lisa, Acta Phys. Polonica B 47, 1847 (2016). 
Jing Yang, Wei-Ning Zhang, Nucl. Sci. Tech. 27 , 147 
(2016). 

R. Ghetti et al., Phys. Rev. Lett. 91, 092701 (2003); G. 
Verde et al., Eur. Phys. J. A, 30, 81 (2006). 

W. G. Gong, W. Bauer, C. K. Gelbke, and S. Pratt, Phys. 
Rev. C 43, 781 (1991). 

Y. G. Ma et al., Phys. Rev. C 73, 014604 (2006). 

N. Colonna et al., Phys. Rev. Lett. 75, 4190 (1995). 

R. Ghetti et al., Phys. Rev. C. 69, 031605(R) (2004). 
L.W. Chen, V. Greco, C. M. Ko, and B. A. Li, Phys. 
Rev. Lett. 90, 162701 (2003). 

M. Lv, Y. G. Ma, G. Q. Zhang, J. H. Chen, D. Q. Fang, 
Phys. Lett. B 733, 105 (2014). 


37 


38 


M. Lv, and Y. G. Ma, J. H. Chen, D. Q. Fang, G. Q. 
Zhang, Phys. Rev. C 95, 024614 (2017). 

B. T. Haar and R. Malfliet, Phys. Rev. C 36,1611(1987); 
G. Q. Li and R. Machleidt, Phys. Rev. C 48, 1702(1993). 
X. Z. Cai, J. Feng, W. Q. Shen, Y. G. Ma et al., Phys. 
Rev. C 58, 572 (1998); J. Y. Liu, W. J. Guo and G. J. 
Wang et al., Phys. Rev. C 86, 975 (2001). 

Y. X. Zhang, Z. X. Li and P. Danielewicz, Phys. Rev. C 
75, 034615 (2007); B. Chen, F. Sammarruca, and C. A. 
Bertulani, Phys. Rev. C 87, 054616 (2013). 

R. Lednicky, Sov. J. Nucl. Phys. 35, 770 (1982). 

R. Lednicky, Phys. Ato. Nucl. 71, 1572 (2008). 

R. Lednicky, Phys. Part. Nucl. 40, 307 (2009). 

S. E. Koonin, Phys. Lett. B 70, 43 (1977). 

R. Lednicky and V.L. Lyuboshitz, Heavy lon Physics 3, 
93 (1996). 

P. Danielewicz, Phys. Lett. 146B, 168 (1984). 

J. Aichelin, A. Rosenhauer, G. Peilert, H. Stoecker, and 
W. Greiner, Phys. Rev. Lett. 58, 1926 (1987). 

T. Z. Yan, Y. G. Ma et al., Phys. Lett. B 638, 50 (2006). 
Y. G. Ma and W. Q. Shen, Phys. Rev. C 51, 710 (1995). 
G. Q. Zhang et al., Phys. Rev. C 84, 034612 (2011); Y. 
Z. Xing et al., Chin. Phys. Lett. 33, 122501 (2016). 

Z. Q. Feng, Nucl. Sci. Tech. 26, 520512 (2015). 

C. L. Zhou et al. , Phys. Rev. C 90, 057601 (2014). 

C. Tao, Y. G. Ma, G. Q. Zhang et al., Nucl. Sci. Tech. 
24, 030502 (2013); Phys. Rev. C 87, 014621 (2013). 


J. Chen, Z. Q. Feng, and J. S. Wang, Nucl. Sci. Tech. 27, 
73 (2016). 

B. S. Huang and Y. G. Ma, Chin. Phys. Lett. 34, 072401 
(2017). 

W. J. Xie, F. S. Zhang, Chin. Phys. Lett. 32, 122502 
(2015). 

T. T. Wang, M. Lv, Y. G. Ma et al., Chin. Phys. Lett. 


32, 062501 (2015). 

S. Pratt, Phys. Rev. Lett. 53, 1219 (1984). 

K. Chen et al., Phys. Rev. 166, 949 (1968). 

E.Fermi, Theor Phys. 5, 570 (1950). 

P. J. Siemens, and J. O. Rasmussen, Phys. Rev. Lett. 42, 
880 (1979). 

E. Schnedermann, J. Sollfrank, and U. Heinz, Phys. Rev. 
C 48, 2462 (1993). 

T. Lefort, D. Doré, D. Cussol, Y.G. Ma, J. Péter, R. 
Dayras et al., Nucl. Phys. A 662, 397 (2000). 
Helgesson,J. et al., Phys. Rev. C 56, 2626-2635 (1997). 
S. Zhang, Y. G. Ma, J. H. Chen, and C. Zhong, Advances 
in High Energy Physics 2016, 9414239 (2016). 


whole-rapidity 


mid-rapidity 


projectile/target rapidity 


" 
1.2} fot, u (a) 1.8 (e) 
0.8} ee, 1.2 | b=3fm 
—a—7=0.2 
—v—n20.5 1 0.6 
—4—=0.7 
eens 0.9 1 0.0 
(b) 1.8 (f) 
1.2 | b=6fm 
1 0.6 
i 1 0.0 
1.2} $ (c) 1.8 (9) 
S eU 
E 0.8 1.2 |b-9tm 
0.41 0.6 
0.0 
1.2} (d) 1.2} (I) 
0.8} 0.8} 1 
b=12fm 
0.41 0.4. 
0.0 ; ; , , , , 10.0 4 
0 20 40 60 80 100 0 20 40 60 380 £100 0 20 40 60 80 100 
q(MeV/c) 


FIG. 2: (Color online) Momentum correlation function of proton-pairs for Au + Au collisions at 1A GeV with different in- 
medium reduction factors. From top panel to bottom panel, each panel corresponds to impact parameters, namely b = 3, 6, 
9 and 12 fm, respectively. From left to right panel, each panel represents the correlation function constructed for the proton 
pairs within whole rapidity, mid-rapidity and projectile/target rapidities, respectively. 


(a) b=3 fm, n=0.2 
red circles: mid-rapidity cut 
blue squares: whole-rapidity 


/2-e—e 


© g 


p 2 
aa Hom 


400 800 1200 1600 
E (A MeV) 


=20MeV/c) 


pa 


PP 
ho 


C 
o 


0 2 4 6 8 10 12 
Impact parameter (fm) 


(c) E= 1A GeV, b=3 fm 


0.0 0.2 0.4 0.6 0.8 1.0 
n 


FIG. 3: (Color online) (a) Cpp strength at 20 MeV/c as a function of the incident energy, which is calculated at an impact 
parameter of b = 3 fm and in-medium reduction factor of 7 = 0.2; (b) Cpp strength at 20 MeV/c as a function of the impact 
parameter, which is calculated at 1A GeV and the in-medium reduction factor 7 = 0.2; (c) Cpp strength at 20 MeV/c as a 
function of the in-medium reduction factor, which is calculated at b = 3 fm and 1A GeV. Note that the red circles corresponds 
to mid-rapidity cut and the blue squares whole-rapidity cut. 
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FIG. 4: (Color online) Gaussian source radius as a function of the incident energy in different rj-factor at fixed impact parameter 
b = 3 fm (left column), as a function of impact parameter with different 7 factors (middle column), as a function of in-medium 
nucleon-nucleon cross section reduction factor at different impact parameter (right column). From upper to lower rows, it 
corresponds to the Gaussian radius of free protons emitted from the whole rapidity window, the mid-rapidity window, and the 
target or projectile rapidity window, respectively. 
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FIG. 5: (Color online) The x-square obtained from the fits of the proton-proton momentum correlation function by the Lednicky 
et al. analytical formalism calculation as a function of the radius of the Gaussian source. Left panels: each panel represents 
different fixed impact parameters with different 7 factor; Right panel: each panels represent different fixed 7 factors at different 
impact parameters, respectively. 
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FIG. 6: (Color online) Left panels: Blast-wave fitting to pr spectra of protons in the mid-rapidity region. upper panel: incident 
energy dependence for b=3 fm and 7=0.2; middle panel: impact parameter dependence at E = 1A GeV and 7 = 0.2; lower 
panel: in-medium cross section dependence at b = 3 fm and 1A GeV. Right panels: the fitted radial flow parameter (3) as 
a function of beam energy (upper panel), impact parameter at 1A GeV and 7 = 0.2 (middle panel), in-medium cross section 
reduction factor (lower panel) at 1A GeV and b = 3 fm. 
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FIG. 7: (Color online) Correlation strength at q = 20 MeV/c as a function of the radial flow parameter (8) in different 
conditions. Upper row: beam energy is a variable as shown in figure, at fixed b = 3 fm and 7 = 0.2 (upper row); Middle 
row: impact parameter is a variable as shown in figure, at fixed E = 1A GeV and 7 = 0.2; Bottom row: the in-medium 
nucleon-nucleon cross section reduction factor is a variable as shown in figure, at fixed E — 1A GeV and b — 3 fm. 


